Astronomy <fe Astrophysics manuscript no. paperbib 


February 2, 2008 


(DOI: will be inserted by hand later) 





Analysis of two-point statistics of cosmic shear: 
II. Optimizing the survey geometry 

Martin Kilbinger^ and Peter Schneider^'^ 

^ Institut f. Astrophysik u. Extr. Forschung, Universitat Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 
^ Max-Planck-Institut f. Astrophysik, Postfach 1317, D-85741 Garching, Germany 

Received / Accepted 

Abstract. We present simulations of a cosmic shear survey and show how the survey geometry influences the 
accuracy of determination of cosmological parameters. We numerically calculate the full covariance matrices Gov 
of two-point statistics of cosmic shear, based on the expressions derived in the first paper of this series. The 
individual terms are compared for two survey geometries with large and small cosmic variance. We use analyses 
based on maximum likelihood of Gov and the Fisher information matrix in order to derive expected constraints 
on cosmological parameters. As an illustrative example, we simulate various survey geometries consisting of 300 
individual fields of 13' x 13' size, placed (semi-)randomly into patches which are assumed to be widely separated on 
the sky and therefore uncorrelated. Using the aperture mass statistics (Mfp^, the optimum survey consists of 10 
patches with 30 images in each patch. If Qm, os and F are supposed to be extracted from the data, the minimum 
variance bounds on these three parameters are 0.17, 0.25 and 0.04 respectively. These variances raise slightly 
when the initial power spectrum index rig is also to be determined from the data. The cosmological constant is 
only poorly constrained. 
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1. Introduction 

' Weak gravitational lensing by the large-scale mat- 
. ter distribution in the Universe, called cosmic shear, 
' has become a valuable tool for cosmology since its 
. first detection in 2000 JBacon et al.ll2000t iKaiser et al.l 
; l2000t Ivan Waerbeke et all l2000t IWittman et all \200(f . 
. Constraints on cosmological parameters, in particular the 
' (dark+luminous) matter density parameter and the 
power spectrum normalization cts, have been obtained 
from cosmic shear wit h survey area s of up to several dozen 
square degrees fe.g. Ijaryis et all l2003t fciadders et all 
2005; 'van Waerbeke et al.ll200li l2002t iMaoli et al.ll20o'ir 
Rcfrcgicr et al...2002t) . 

Cosmic shear probes the statistical properties of the 
total matter distribution projected along the line-of-siglit 
between the observer and distant galaxies which are typ- 
ically at redshifts between 0.5 and 2. It is independent 
of any possible bias between dark and luminous matter 
— e.g. no assumptions about how galaxies trace the dark 
matter have to be made. 

Cosmic shear is sensitive to a large number of cos- 
mological parameters, most notably on Q^^as and the 
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shape parameter F, but also to the source galaxy red- 
shift distribution. The dcpendancy on these parameters 
is partially degenerate. Some of these near-degeneracies 
can be broken when weak lensing is combined with other 
cosmological measurements like CMB anisotropy experi- 
ments, the statistics of the Lyman-a forest or galaxy red- 
shift surveys. The parameter dependencies are very dif- 
ferent for the individual methods, for example the fl^-as- 
degene racy is nearly orthogonal b etween cosmic shear and 
CMB llvan Waerbeke et"al]l2002() . Even the most precise 
measurement of cosmological parameters u p to now, which 
comes from the first data release of WMAP ijBennett et alJ 
l200.lHSDergel et al.ll200.ll) . can be i mproved substantially 
when wea k Icnsingdata is added ijContaldi et al.l 120031: 
IHu fc Tee mark 19^. 

In the first paper of this series ijSchneider et al.l|2002L 
hereafter Paper I), we reviewed the properties and rela- 
tions of various two-point statistics of cosmic shear. We 
defined unbiased estimators and calculated their covari- 
ances. In Sect. 5 of Paper I, an explicit expression for the 
covariance of the two-point correlation function was given. 
Using that, we calculated expected constraints on various 
cosmological parameters for a cosmic shear survey. 

However, this ansatz is an approximation which is only 
valid for a large and connected survey area. Any real cos- 
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mic shear survey will most likely consist of single uncon- 
nected fields-of-view and have a complicated geometry. In 
this paper, we present a method which allows one to cal- 
culate the covariance of the two-point functions of cosmic 
shear for an arbitrary survey geometry. 

The measurement of cosmic shear with a sufficient 
high precision to constrain cosmological parameters re- 
quires many independent lines-of-sight, lowering the sam- 
pling variance ("cosmic variance"). On the other hand, 
it is important to measure the shear on a large range 
of angular scales. Even with modern wide-field imaging 
cameras, separations of more than a few degrees cannot 
be accessed by individual fields-of-view — one has to ob- 
serve some fields near to each other and measure galaxy 
shape correlations across individual fields. In a recent work 
()jarvis et alJl2003|) . cosmic shear has been measured from 
1 to 100 arc minutes, and constraints on f2m and cts have 
b e en obtained. 

iKaiserl (^2^ remarked that for a survey of 9 square 
degrees, consisting of a single 3° x 3°-field. the noise due to 
the intrinsic cUipticity dispersion can be neglected on large 
scales, because of the huge number of galaxy pairs. He also 
noted that cosmic variance is dominant, and that "sparse 
sampling", meaning the distribution of smaller fields on 
a larger region of the sky, reduces the cosmic variance 
dramatically. 

A cosmic shear survey has to cover a large area con- 
taining hundreds of thousands of galaxies, whose shapes 
can be determined. Because telescope time is limited, one 
has to carefully choose the locations of the pointings, in 
other words, the geometry of the survey. In this work, 
cosmic shear surveys with different geometrical configura- 
tions are simulated. They are compared with respect to 
their ability to constrain cosmological parameters, using 
two-point statistics. 

The individual fields-of-view of the simulated surveys 
are placed in patches on the sky in order to measure the 
shear correlation on large angular scales. Several patches, 
distributed randomly on the sky in order to reduce cosmic 
variance, build up the survey. 

The comparison of the geometries is done with a likeli- 
hood analysis using the covarianccs of two-point statistics 
of cosmic shear. These were derived in Paper I, Sects. 4 
and 6. In this paper, the covariance matrices are calcu- 
lated via a Monte-Carlo-like method using the simulated 
galaxy positions of the surveys. 

The different survey geometries considered here are 
presented in Sect. El In Sect. |3 we review the two-point 
statistics of cosmic shear relevant for this work, as well 
as their estimators and covarianccs, as derived in Paper 
I. The method for the numerical calculation of the covari- 
anccs is given in Sect. 0] Results for some patch geome- 
tries are presented in Sect. 14.31 A likelihood analysis is 
performed in Sect. [S] where the expected constraints on 
pairs of cosmological parameters is considered and com- 
pared for various cosmic shear survey geometries. Finally, 
in Sec. El the Fisher information matrix is used to compare 
constraints on three and four parameters simultanously. 



This work is intended to be a preparation for a cosmic 
shear survey with the wide-field camera VIMOS on the 
8.2m ESO-VLT telescope UT3 (Melipal). The numerical 
codes used for this paper are publically available^. 

2. Survey geometries 

The different survey geometries considered in this work 
consist of circular patches, in which individual images 
are distributed randomly, but non-overlapping. We define 
"image" as one single field-of-view. It is assumed that the 
shear correlation functions can be measured across image 
boundaries, thus the cosmic shear can be determined in 
principle on scales up to the patch diameter. Because there 
are always bright stars or foreground galaxies which have 
to be avoided, we cannot specify in much more detail the 
image positions, thus a random distribution of the images 
in a patch is assumed. 

The patches are assumed to be separated by at least 
several degrees. On these scales, the correlation functions 
are virually zero, so different patches can be considered as 
uncorrelated; they probe statistically independent parts of 
the large-scale structure. 

In our simulations, the individual images are 13' x 
13'-fields, corresponding roughly to the field-of-view of 
VIMOS. A survey consists of P patches of radius R, each 
patch containing N images. The total number of images, 
n = P ■ N is kept fixed for all geometries. The larger the 
number of patches P, the smaller is the cosmic variance. 
On the other hand, the larger the number of images N 
per patch, the larger is the number of galaxy pairs for 
which the correlation is measured, thus the lower is the 
shot noise. One of the goals of this work is to find a con- 
figuration which is optimal in the sense that the two-point 
correlation function can be measured most accurately; we 
characterize this 'accuracy' by considering constraints on 
pairs of cosmological parameters from the measurement of 
the correlation function. 

We use a total image number of n = 300, correspond- 
ing to 14 square degrees for the whole survey. For N, be- 
ing a factor of n, the values 10, 20, 30, 50, and 60 are 
considered, corresponding to geometries with P ~ 30, 15, 
10, 6 and 5 patches, respectively. An illustration of some 
patches is given in Fig. ^ 

The patch geometries are compared to a configuration 
which consists of 300 single uncorrelated images, where 
'uncorrelated' again means separated by at least several 
degrees. This configuration has the smallest possible cos- 
mic variance, but the shear correlation can be only mea- 
sured up to y/2 ■ 13 arc minutes. 

3. Two-point statistics of cosmic shear, their 
estimators and covariances 

We fo llow here the notation of iBartelmann fc Schneide3 
(|200l|) . See Paper I for a more detailed presentation of 

^ http : //www . astro . uni-bonn . de/~kilbinge/ cosmicshear 
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Fig. 1. Realisations of patches; the squares represent the 
individual images. N is the number of images per patch, R 
the patch radius. Each image has a size of 13 arc minutes. 

the formulae reviewed in this Section. Two points i9 and 
"d + 6 define in a natural way a direction given by the 
connecting vector 9, with respect to which the tangential 
and cross-component of the shear, 



7t = - 3? (76" 



-2i.p\ 



and 



7> 



= - 3 (76-2*'^) 



(1) 



are defined, where (p is the polar angle of the connecting 
vector. The two independent shear correlation functions 
C±(i^) are related to the power spectrum of the projected 
matter density P„, 



(7tW7t(i9 + 0))±(7xW7x(i5 + e)> 
d«P«(£)Jo,4(^e), 



1 

2^ 



(2) 
(3) 



where the first-kind Bessel function Jo (J4) corresponds t o 
the ('— ') correlation function, see e.g. iKaisCT ljl992t) . 
Another second-order statistics of cosmic shear is the dis- 
persion of the weig hted tangential s hear in a circular aper- 
ture B of radius 9 (ISchneidelEiai) , 



This aperture mass measures the gradient component or 
E-mode of the shear. Analogously to Q), one can also 
define the weighted cross-component of the shear in an 
aperture, {M±), which is a measure of the curl component 
or B-mode only. 

The two aperture mass statistics can be expressed in 
terms of the two-point correlation functions. 



28 



2612 



de' e' 



±^_(6i')r_ 



(6) 



where and T_ are given explicitly in ISchneider et alJ 
(l2002t) . 

The shear due to the tidal gravitational field of the 
large -scale structure is a pure gradient field^ ijKaiserl 
therefore, no B-modes should be present. However, 
there are other effects which can produce B-modes; these 
are mainly systematic measurement errors and intrinsic 
galaxy orientation correlations. The contribution from the 
latter can be reduced if the survey is deep, or if photomet- 
ric redshi ft information of th e sourc e galaxies is taken into 
accou nt I King fc Schneider! I2OO2I: iHevmans fc HeavensI 
l200.3h . Minor contribution s to B-modes are source cluster- 
ing llSchneider et al.l2002l) and higher-order lensing effects 
l)jain et al.ll2000|) . 

The separate measurement of E- and B-modes allows 
one to quantify contaminations to the gravitational shear 
signal. In recent cosmic shear surveys, a non-zero B-mode 
signal has been rneasured llvan Waerbeke et~al] l2002t 
iHoekstra et al.ll2002t I.Tarvis et al.ll200.1^■ using the statis- 
tics (jnj. However, a recent re-analysis of the VIRMOS- 
DESCART data shows no significant B-mode signal any 
more; the previously found B-modes were obviously due 
to insufficient PSF corrections'^. 

3.1. Shear estimators 

The shear estimators used here are similar to those intro- 
duced in Paper I. For simplicity, all weight factors which 
account for differences in the precision of ellipticity mea- 
surement of individual galaxies are set to unity. Further, 
we assume that the correlation function is to be estimated 
in logarithmic bins; therefore, the following function is de- 
fined, 



Map(0)= j d^i?g(z?)7tW. 

B{e) 



(4) 



With the weight function Q{§) = 2/^^ dp pU{p) - 
U{'d), where U is a compensated filter function, i.e. 
J^di}dU{i3) = l|Schneider et al.lll998|) . its dispersion 
is related to the power spectrum by 



diep^ie) 



m' J 



(5) 



for |lni9-ln6l| < 
otherwise 



(7) 



which selects bins with logarithmic bin width A In 6 
around z?. Then, if the i-th galaxy is located at angular 
position 9i with the observed ellipticity e,;, 

l±W = ^^^(e.te.t ±£^xe,x)A4|e. - ^,1), (8) 



^ at least in first approximation 
van Waerbeke, priv. comm. 
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are unbiased estimators of the correlation functions 
Here, A^p {'&) is the number of galaxy pairs in the bin cor- 
responding to I?. The double sum is performed over all 
galaxy pairs. 

Unbiased estimators for the aperture mass dispersions 
(A/4) and (Mj) are 



M±{9) 



26*2 



1=1 



,(9) 



where an index has been attached to the bin width A'&i 
to account for variable (e.g. logarithmic) bin widths. The 
limit / of the sum must be chosen such that 9 is half the 
upper limit of the I-th bin. 

3.2. Covariances 

In Paper I we calculated the covariance matrices of the 
estimators defined in the last section. These consist of 
several terms which we call shot noise or diagonal term 
(Z?), mixed term (A/) and pure cosmic variance term {V). 
We perform the following decomposition: 

Cov(e+,7?i;e+,^2) = D + M++ + V++ 
Cov(^"_,?9i;|_,?92) = D + + (10) 

where the individual terms are 



r 

ijk 



V++ = 



|^A^^A^^C+(zOe+(jfc) 

ijkl 

^ ^ ^l^^ cos4((^,, - ^,,)^+{3k) (11) 



ijk 



= ^ E ^^^i^^t ( cos4((/?y - ipa - ipjk + fki) 

ijkl 

X + cos4(^y - (Pki)^+{il)^+{jk) 

^ ^ Y^AlAl cos4(^,, - ^jk)^Ajk) 



ijk 



F 

ijkl 

with F = iVp(i9i)iVp(z?2), A^^' = A^i\e,~9,\) and 
= ^±(|f» - for brevity S{^i - ^2) is unity if 
"i?! and 1^2 are in the same bin and zero otherwise, is 
the ellipticity dispersion of the galaxies in the absence of 
shear. 



The covariances of the aperture mass dispersions are 
Cov(A<±;0i,02) = T^EE^'^'^^^^'^^- 



-1- ^ i=l j=l 



T. 



Cov(^+,z?,;e+,^,) 



r-(|)T_('^)Cov(^_,^.;C-,^,) 



±r+(^')T_('^)Cov(e+,^,;e-,i?0 



(12) 



The upper limits /fc, k = 1,2, must be chosen such that 6k 
is half the upper limit of the /fe-th bin. For brevity, the no- 
tions Cov+_|- = Cov(^+, .; .), Cov__ = Cov(^_, .; .) 

and Cov^ = Cov(^+, .; .) are used from now on. 

We note here that the expressions for the cosmic vari- 
ance terms V of the covariances are only valid if the 
shear field is Gaussian. On scales below ~ 10 arc min- 
utes, the non-Gaus sianity of the shear field ge ts impor- 
tant, e.g. Fig. 4 of Ivan Waerbeke et al.l l)2002() . see also 
IScoccimarro et all l)l999|) . On scales below 1 arc minute, 
the shot noise term D dominates over V, see Fig. 0] of 
this paper and Fig. 3 of Paper I, thus with the Gaussian 
assumption we expect to slightly underestimate the co- 
variances in this transition regime between 1 and 10 arc 
minutes. 



4. Numerical calculation of the covariance 
matrices 

Given a model for the shear correlation functions or, 
equivalently, the power spectrum P^, the covariance ma- 
trices (|10ll only depend on the positions of the observed 
galaxies, in other words, on the survey geometry. For a 
given data set with known positions of observed galax- 
ies, it is straightforward to calculate the covariances. An 
a priori estimate of the covariances is made using simu- 
lated galaxy positions for the summations in (|11() . Note 
that only the positions of the galaxies have to be sim- 
ulated, not their ellipticities. In order not to introduce 
artificial Poisson noise, the galaxies are not distributed 
ra ndomly but subran domly onto the fields, see Chapter 7 
of lPress et all (Il992l) . 

Throughout, 20 logarithmic bins in angular separation 
are used, the smallest bin being centred around 10 arc 
seconds. The largest separation considered is cither ^/2-13' 
for the uncorrclated images, or equal to the patch radius 
R in the case of a patch geometry. Thus, the bin widths 
differ for geometries with different patch radii. 

The number of galaxy pairs per bin is shown in Fig. 
121 for a single 13' x 13' -field. For intermediate angles, 
Np{'d) (X '0'^, as follows from eq. (26) of Paper I, with 
A-d (X 'd for logarithmic bins. Deviations show up for large 
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Fig. 2. Number of galaxy pairs per angular bin A'p against 
the bin centre for a single 13' x 13' -field. The points 
are the pair numbers obtained from the subrandom dis- 
tribution, the dashed line is the approximation A^(i?) = 
l-nAr^'d^'d (Paper I), where A is the survey area and 
n the galaxy number density, which was taken to be 
30 per square arc min ute. The dotted curve is from 
iKins- SchneideH 



the triple sums and 1 000 pairs per bin for the quadruple 
sums. 

The calculation of the individual addends in the co- 
variances is straightforward. The correlation functions are 
obtained by linear interpolation between grid points cal- 
culated beforehand, using (j2Jl. The cosine of the sum of an- 
gles is expanded into a sum over products of cosines, and 
by using the relation cos4(^ = 2^1x^1 (x^ + y'^Y ^ 1]^ ~ Ij 
where is the polar angle of the vector {x,y), no single 
time-consuming trigonometric function actually has to be 
evaluated. 

4.1. Cosmological model 

For the power spectrum of the matter fluctuations, we 
assume an initial power spectrum Pi cx fc"°, the trans- 
fer fu nction for Cold Dark Matter from iBardeen et alJ 
and the fitting formula for the non-linear evolution 
of lPeacock fc Doddd lll996tl. The redsh ift distribution of 
the source galaxies is llSmail et al.lll99,ql) . 



p{z)dz 



/3 



- e-(-/-°)'d2 



(13) 



zor(3//?) v^o. 
where F denotes the Eulerian gamma function. Our ref- 
erence cosmology is a flat ACDM model with ft^ ~ 0.3, 
f^A = 0.7, the shape parameter F = 0.21, rig = 1 and the 
normalisation erg = 1. The parameters of the redshift dis- 
tribution are zq = 1 and /? = 1.5, which corresponds to a 
mean source redshift of « 1.5. 



scales, which are due to boundary effects and for very 
small scales, because of the subrandom galaxy distribu- 
tion. The curve agrees well w ith the theoretical expecta- 
tion ((King & Schneiderll2n03() . 

Throughout, we set the intrinsic ellipticity dispersion 
to (7j = 0.3. The number of galaxies per square arc minute, 
for which a shape mesurement is feasible, is set to 30. To 
obtain this number density, a limiting i?-band magnitude 
of about 25.5 is needed. 

Because the number of operations for the calculation 
of the covariances increases with the number of galaxies to 
the fourth power, it is not feasible, except for a very small 
survey area, to sum over all galaxy positions. Instead, a 
random subsample of galaxies is used. 

The summations over galaxy positions can be written 
as sums over pairs of galaxies, which have a flxed separa- 
tion for each matrix element (t?!, •&2), as determined by the 
A- functions. Thus it is convenient to store galaxy pairs for 
each angular bin. The mixed terms M (|ll(l can then be 
split up into a sum over all those §1- and 'd2-j>SLVcs which 
have a galaxy in common. This can be done efficiently if 
the pairs for each bin are sorted by galaxy number. The 
cosmic variance terms V are simply double sums over all 
t?i- and i92-pairs. 

In order to decrease the computing time to a feasible 
value, we use random subsamples of all galaxy pairs for 
the summations, which consists of 300 000 pairs per bin for 



4.2. The covariance matrices 

In Fig. 13 contour plots of the three covariance matri- 
ces arc shown for a single 13' x 13' -field. The diagonal 
of Cov_|_+ is enhanced above the non-diagonal elements 
only for small angular separations, which is mainly due 
to the (l/A'p)-dependence of the shot noise term D. The 

Gov -matrix shows a more contrasted diagonal, with a 

rapid falloff away from the diagonal. This stems from the 
fact that ^_ filters the convergence power spectrum more 
locally than ^-|_, resulting in smaller intercorrelations of 

different angular scales. The asymmetric matrix Cov^ 

has also negative diagonal elements, meaning anticorrela- 

tion. For ~ 1-4 Cov^ is zero. 

The covariance matrices of the two aperture mass 
statistics arc quite similar to each other, indicating that 
the third and fourth terms in (|12|l are small in compar- 
ison to the first two terms. Moreover, they also resem- 
ble Gov , which is due to the resembling functional be- 
haviour of and (A/^p); the filter fimctions for both 
statistics (^4(2;) and [J4{x)/x]'^ respectively) are strongly 
peaked and thus filter the power spectrum very locally. 

4.3. Comparison of various geometries 

In Fig. 01 the individual terms (|ll|l contributing to the 
diagonal (t9i = 7^2) of the covariance matrices (IIOII . calcu- 
lated for four survey geometries, are compared. One sees 
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Fig. 3. Contour plots of the covariancc matrices of the correlation functions and the aperture mass dispersions 
(O (Mfp) ("E-mode") and (M^) ("B-mode"), for a single 13' x 13' -image. Logarithmic contours are used, ranging 
from 2 X 10"^"^ to 10~^, with 4 spacings per decade. Negative values are plotted as dotted contours. The 'boxiness' 
of some contour lines is due to the coarse resolution; the matrices are calculated on a logarithmic 20 x 20 -grid for 
and '(?2- The axes show iJi and in arc minutes. Note that the aperture mass covariances can only be calculated up 
to half the largest scale where the covariances of the correlation functions are known. 



the (l/Afp)-dependence of D; it decreases as the number of 
pairs increases towards larger separations. Only for scales 
comparable to the image or patch boundary, the number 
of pairs decreases. This leads to an increase of D (see also 
Fig- El- For Cov+-)_, the cosmic variance V dominates on 
scales larger than about one arc minute. Cov is domi- 
nated by D. The Cov-| -elements are very small compared 

to those of the two symmetric matrices. 

In the upper row of Fig.01 all terms for the two extreme 
geometries regarding cosmic variance are shown: the 300 
uncorrelated images, and the configurations with only five 
patches with a radius of 80 arc minutes. The difference is 
about a factor of two in the cosmic variance term V+-)- 
The other terms arc quite similar for the two geometries, 

except on scales comparable to the image size. V is 

much smaller; at this level, the differences between the 
geometries are presumably mainly due to numerical noise, 
as well as the negative value at about 0.6 arc minutes for 
the patch geometry . Thus, ^+ is much more affected by 
cosmic variance than ^_ . 

Obviously, the covariances corresponding to the patch 
geometry extend to larger angular scales than those for 
the uncorrelated images. One important question which is 
adressed in this paper is whether the additional informa- 
tion of the shear on large scales can compensate for the 
larger cosmic variance on smaller scales. 



The lower row of panels of Fig. 21 compares the co- 
variance terms of patch geometries with the same radius, 
but with a low (TV = 10) and a high {N = 60) image 
density in the patches. In the first case, there is quite a 
sharp transition at a scale where the image boundary is 
exceeded, all terms increase at about 10 arc minutes. The 
case iV = 60 shows a less drastic change; because of the 
higher image density, a large number of pairs at this sep- 
aration on different images is found, thus the transition is 
smoother. 

4.4. Comparison with the approximation of Paper I 

We compare the single terms of the covariance matrix 
Cov4_-(- H10(l obtained by the summation of simulated 
galaxy positions presented in this paper with the inte- 
gration method from Paper I, Sect. 5, for a large con- 
nected field. We found deviations at large separations, 
where the assumption validating the approximation for 
the integration method breaks down. Further, for small 
separations, the diagonal elements of the mixed term M^^ 
term are enhanced in the summation method. This is due 
to the discreteness of the galaxies: For di = i?2, there are 
Np{'di) summands where j = k. This gives a contribution 
of iVp(i?i) ^+(0) which is not present for the off-diagonal 
elements. A similar but smaller contribution adds to V++. 
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Fig. 4. The diagonal elements of the covariancc matrices (|10|l . split up into the individual terms D, M and V 
Upper row: 300 single uncorrelated images, where the largest scale is \/2 ■ 13' (thick lines) and a patch geometry with 
iV = 60 and = 80' (thin lines). Lower row: N = 10, R = 100' (thick Unes) and N = 60., R ^ 100' (thin lines). 
The left and middle panels are logarithmic plots, the axis of the right panel is linear. Note that the bin widths 
are different for the different geometries which affects D. Because logarithmic bins are used, the differences in the bin 
widths are small, however. 



These effects do not enter the integration method, where 
a smooth galaxy distribution is assumed. 

Further differences between the two methods were 
found for the cosine-part of the V+-|--term which is due to 
fact that f _ decreases very slowly for large separations and 
thus the (/3-integration in eq. (34) of Paper I obtains a con- 
siderable contribution from separations which are larger 
than the field boundary. 

Altogether, the deviations are quite small. In particu- 
lar, the resulting likelihood contour plots (see next section) 
are very similar, which confirms the consistency of the two 
methods. 



5. Likelihood analysis 

By using the covariance matrices, we construct an a pri- 
ori estimate on cosmological parameter constraints from 
a cosmic shear survey. This allows us to compare differ- 



ent survey geometries. As in Paper I, we use the following 
figure-of- merit : 



where the superscript ' denotes the fiducial model and 
p is a set of cosmological parameters which is tested 
against the fiducial model, as specified in Sect. 14.11 The 
summation indices label the angular bins of the corre- 
lation functions. As noted in Paper I, either ^+ or ^_ 
can be inserted for ^ in (|14|l . in which case the corre- 
sponding covariance matrices Cov_)_4_ or Gov , respec- 
tively, have to be used for Gov. The resulting function 
is called x+ or X- respectively. The figure-of-merit Xtot 
combining all information is obtained using the vector 
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^ = (C+i: ■ • ■ , S.+N, £.-n) together with the 

block matrix 



Gov 



Cov+_|_ Cov_|-_ 



Cov^_ Cov_ 



(15) 



Analogously to l(TH) . a figure-of-merit using the (Mfp)- 
statistics can be defined, 



X {{Ml{e,)) - {Ml{e,))') , (16) 

which is a likelihood measure using only the E- modes. 
Unfortunately, Xsi using M± and thus testing the B- 
modes, cannot properly be defined in this way. In the 
model used throughout in this work, no B-modes are 
present, thus {M^} = 0. The correlation functions Q 
only depend on the E-mode power spectrum; up to now, 
no models for the B-mode power spectrum have been ob- 
tained. Predictions of the amplitude of intrinsic galaxy 
alignment differ by orders of magnitudes and although 
some observational results have been presen ted, stringent 
const raints have not yet been obtained fsee iBrown et all 
I2OO2I and references therein). 

Two of the four parameters r2ni,(T8,r and zq a-re var- 
ied in a few combinations, while all others are kept fixed, 
with the exception that a flat universe is assumed through- 
out. In Figs. [SI and ini the expected likelihood contours for 
the two extreme geometrical configurations regarding cos- 
mic variance are plotted: the uncorrelated images and the 
{N = 60, i? = 80')-patch geometry. 

5.1. and dg 

The strongest degeneracy between two parameters exists 
for f2ni and erg, which is also expect ed from simulations 
and c osmic shear measurements fe.g. lvan Waerbeke et al.l 
I2OO2I) . In order to compensate for the high elongation of 
the contours, we use the combined parameter 



= (78 



0.41 + 0.59 



03 



-0.68 



(17) 



which has been obtained in Paper I by fitting the minimum 
"valley" in the (rini-o'8)-plot. Clearly, the x?--contours are 
more extended in the case of the uncorrelated images than 
for the patch geometry. This is because ^_ contains much 
information on large scales which is absent in the case of 
the uncorrelated images. In contrast, the x+-contours are 
tighter in this case than for the patch geometry. 

Furthermore, in both cases, the difference between 
and x^oj are small. Thus, most of the information concern- 
ing cosmology is contained in the additional informa- 
tion coming from ^_ is relatively small. 



5.2. Other combinations 

There is also a strong degeneracy between other combina- 
tions of parameters, as seen in Figs. [S] and In all cases, 
the Xtot" a-iid x+-contours are tighter for the uncorrelated 
images, whereas the opposite is true for X-- The patch 
survey geometry yields constraints from ^_ compatible to 
those from ^_|_ , in particular when zq is one of the parame- 
ters; for the combination CTf^-ZQ^ the x--contours are even 
tighter than the X-i--contours. 



5.3. Quadrupole moments 

For a more detailed analysis, a quantitative description 
of the x^-contour plots presented in the last section is 
needed. The quadrupole moments of the underlying prob- 
ability function can be used as a measure of the surface of 
the contours. These are defined as 



^ /d^pL(pi,p2) jPt-pl) {Pj ~p)) 



(18) 



for i,j — 1,2. The integration variables pi,P2 are the pa- 
rameters which are plotted on the axis in Figs. |31 and El 
pI denotes the value of the parameter of the fiducial set, 
L (X exp(— x^/2) is the likelihood function. 

The determinant of the quadrupole moment is then a 
measure of the surface enclosed by the likelihood contours. 



\/ det Qi 



Q11Q22 — Qr 



12- 



(19) 



If the x^-function was quadratic, the likelihood L would 
be a bivariate Gaussian. If in addition the two parameters 
Pi and p2 were uncorrelated, q oc Api ■ Ap2, where Api is 
the variance of pi. In general, the smaller the value of q, 
the tighter are the likelihood contours, and the better the 
cosmological parameters can be constrained. For the like- 
lihood function L in H18(l . or equivalently, for x^, the three 
functions x+iX- and Xtot '^^^ inserted. The resulting 
g's will be refered to as (74-, q^ and qtot respectively. When 
(|16|l is used as the figure-of-merit, the inferred contour 
surface measure is q^- 

Values of q for different geometries arc given in Fig.^l 
The likelihood contours of all three combinations of pa- 
rameters considered here show the same behaviour. Note 
that the values of q only have a sensible meaning when 
compared to each other for the same combination of cos- 
mological parameters. 

As expected, the tightest constraints are obtained for 
gtot, as it combines the information of all measurements, 
gtot is almost monotonically decreasing with increasing 
patch radius R, and with decreasing number of images per 
patch N. This indicates that cosmic variance is the most 
crucial source of errors: therefore, large patches and most 
notably a large number of directions on the sky should be 
observed. 

q^ behaves similarly to qtot, but gives less tight con- 
straints on the parameters, as was already seen in the con- 
tour plots. The (7_ -behaviour is opposite to gtot and g+. 
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Fig. 5. l-cr, 2-cr and confidence contours of the figurc-of-merit (|14|) . Solid lines correspond to Xt 
xX and dashed lines to The left panels show the results for the 300 uncorrelated images, the right panels are for 
a {N = 60, R = 80')-patch geometry. The parameter Ss « as is defined in H17(l . 



The fact that g_ increases towards small image numbers 
N per patch is due to the weak dependence on cosmic vari- 
ance for Gov Further, small patches give better con- 
straints, indicating that dense sampling on medium scales 
is more important than large-scale information. 

The quadrupole moment measure qe of the (M^p)- 
statistics shows a behaviour similar to g_ . This is because 



of the similarity between ^_ and (A/fp). However, while 
the i?-dependence of is quite monotonic, this is not the 
case for N. For most patch radii R, the tightest constraints 
are archieved for TV = 30. 

For the uncorrelated images, values of q are given in 
Table ^ gtot is smaller than for every patch geometry. 
As expected, g_ is larger than for any patch geometry. 
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Fig. 6. Confidence levels of the figure-of- merit lfn|) . As in Fig.O l-a, 2-a and 2>-a contours are plotted; the solid hnes 
display , the dotted and dashed lines correspond to A- 1 respectively. The left panels show the results for 

the 300 uncorrelated images, the right panels are for a [N = 60, R = 80')-patch geometry. 



Table 1. Quantitative measures of the Vt^-ag, likelihood 
contours for the 300 uncorrelated 13' x 13' images. Table 
entries are in units of 10^^. 



Qtot 


q+ 


q- 




3.825 


4.60e 


26.53 


12.31 



at the minimum of x^. Its eigenvalues and (pairwise or- 
thogonal) eigenvectors can be interpreted as the axes of 
an ellipsoid which determines how fast the log-likelihood 
falls off the maximum in different directions. According 
to the Cramer-Rao inequality, the smallest possible vari- 
ance for any unbiased estimator of a parameter pi, if all 
parameters are to be estimated from the data, is 



because is important on medium scales which arc not 
sampled by the uncorrelated images. 

In Sect. 16.21 we take into account the simultaneous 
determination of three and four cosmological parameters. 



6. Fisher matrix analysis 

The Fisher information matrix is defined as 
a2[-lnL]\ 1 / d'^x^ 



dptdpj 



2 \ dpidpj 



where L is the likelihood function and is the figure-of- 
merit (I14II which depends on the parameter s pi (see e.g. 
iKendall fc Stuartlll969trTegmark et allllOOTl) . The inverse 
of the Fisher matrix is a local measure of the curvature 



(21) 



thus defining a minimum variance bound. 

Because we will later use the (Mfp)-statistics, we in- 
sert (I16II into 120|) and evaluate the equation at the mini- 
mum of to get 



(20) F,, =Y,[Goy-\M + )]H- 



^{Ml,),^{Ml;)^ 



kl 



dpi 



dpj 



(22) 



Given the covariance matrices, we can easily calculate the 
minimum variance for the cosmological parameters used 
in this analysis. 
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Fig. 7. The parameter q H19() for the fim-Eg-contours. On 
the x-axcs, the patch radius R in arc minutes is plotted, 
thus each point represents a patch geometry with N im- 
ages in patches of radius R. The horizontal line in the 
upper panels indicates the value for the 300 uncorrelated 
images (omitted for q- and qe, see also Table 0. 



6.1. Comparison with the likelihood 

First, we compare the minimum variance bound 121|) with 
the likelihood contours from Sect. [Sj using Xtot- this 
case, the Fisher information matrix reduces to a 2 x 2 
matrix. This comparison is shown in Fig. |H1 As expected 
from the Cramer-Rao inequality, the likelihood contours 
are larger than the 1-cr-ellipse from the Fisher matrix. The 
orientation of the Fisher error ellipse coincides with the 
likelihood shapes, i.e. the direction of the minimal and 
maximal degeneracy of parameters is recovered. The larger 
the degeneracy between two parameters, the larger is the 
deviation between the local approximation by the Fisher 
matrix and the likelihood function. For the case where 
the curvature is fixed (flat Universe, ri,n + = 1), the 
degeneracy is much larger than for a fixed cosmological 
constant Ha = 0.7. 



Fig. 8. l-(T-likelihood contours (solid lines) using xtot l|14|l 
compared with the l-a-error ellipse from the Fisher matrix 
(|22|l : the dashed ellipse is for a fiat Universe (as it is the 
case for the likelihood contours), the dotted one is for 
r^A = 0.7. a/h is the axis ratio of the ellipses (the case 
= 0.7 is in parentheses). The configuration is a survey 
of 300 uncorrelated images. 



6.2. More parameters 

Next, we calculate the minimum variance bound of three 
and more parameters out of (fini, fg, F, f^A, "-s) simultane- 
ously, corresponding to a full marginalization over these 
parameters. As seen in Fig. [7| the aperture mass clearly 
gives less tight constraints than the combined two-point 
correlation functions. However, this difference gets smaller 
the more parameters are included, in some cases (for F 
and rig), the minimum variance bound is even smaller for 
(Mfp). Another advantage of the aperture mass is its abil- 
ity to separate E- from B-modes (see Sect.|31l. Therefore, 
we will focus on this statistics from now on. However, we 
must note that because of the local filtering of the power 
spectrum, very large scales are not well sampled by (Mfp). 

Figs. 151- llll show the minimum variance for a different 
number of free cosmological parameters for the individual 
patch geometries. The fixed parameters are set to the val- 
ues given in Sect. 14. ll In the cases where VIa is not a free 
parameter, the prior is a flat Universe (f^A = 1 — ^^m)- 

When taking into account three or more parame- 
ters, the uncorrelated image configuration give very poor 
constraints on these parameters. The minimum variance 
bound is in most cases more than double the value of the 
least optimal patch geometry. The reason is that the lack 
of large-scale information highly raises the degeneracy be- 
tween parameters. This cannot be compensated by the 
small cosmic variance. 

When adding Q.a as a free parameter (compare Fig. 
iniwith Fig. 110(1 . the variance of the shape parameter in- 
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Table 2. The minimum variance for several combinations 
of parameters, for the 300 uncorrelated image configura- 
tion. In each row, those parameters which have an entry 
are assumed to be determined from the data, the other 
parameters are fixed. The counterpart of the three rows 
for the patch geometries are the Figs. 1^-1111 
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r 




0.53 


0.76 


0.09 




0.53 


0.76 


0.14 


0.64 


0.53 


0.77 


0.16 


0.20 



creases by more than a factor of two, whereas the variances 
of f2m and CTg are only slightly enhanced. For large patches, 
the variances of the three parameters are just rescaled, 
whereas for small patches, the change is more complicated. 
The constraint on the cosmolog ical constant is very poor , 
confirming a statement made bv lBernardeau et al 

When the spectral index rig is added (Fig. [TT|l . the 
minimum variance bound of and tJg again increase in 
the same way as when adding SIa, although by a greater 
amount. The variance on F changes completely, taking a 
similar functional form on R and N as the variance on Us . 
The reason for this is, that both parameters determine 
more or less the shape of the power spectrum, whereas 
Qrn and erg influence mainly its amplitude. 

For each data point in Figs. El -El corresponding to 
a survey with N images in P = 300 /N patches, only one 
realization of the random image positions for each of the P 
patches was used. We produced some more realizations for 
two of the patch geometries and found that the scattering 
of the minimum variance bound is about one percent. 

From Fig. (51 the best geometry is a survey with five 
(TV = 60,R= 120')-patches. Considering Fig.lHl a config- 
uration with A'^ = 30 and small R yields the best minimum 
variance bounds. For both cases, a survey with 30 images 
and a patch radius of around 100 arc minutes seems to be 
a good choice. However, the patch radius only has a small 
influence on the minimum variance bound, more impor- 
tant are the number of images per patch and therefore 
the number of patches. 

The difference in the minimum variance bound be- 
tween individual patch geometries can make up to 25 per- 
cent. 



7. Conclusions 

We calculated numerically the covariance matrices (1111 - 
I12() for the second-order estimators of cosmic shear ^± and 
(M^p), which were derived in Paper I, via a Monte-Carlo- 
like technique. Galaxy positions were simulated for various 
cosmic shear survey geometries of 14 square degree area. 
These surveys consisted of a total of 300 images of size 
13' X 13' which were randomly distributed in patches on 
the sky. A number of (semi-)random patch configurations 
were compared to a survey consisting of 300 completely 
uncorrelated images. We performed several analyses based 
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Fig. 9. The minimum variance for the parameters r^m, T 
and (Tg using (j^U and (j^U, for the (A/fp)-statistics. ft^, F 
and (Tg are assumed to be determined from the data, all 
other parameters are kept fixed and a flat Universe is as- 
sumed. Each point in the plot represents a patch geometry 
with N images in patches of radius R. The minimum vari- 
ance bound for the uncorrelated images is given in Table 

HI 
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Fig. 10. The minimum variance for the parameters 
rini,r, (Tg and f^A using the (Mfp)-statistics. See Fig. j^l 
for more details. The minimum variance bound for the 
uncorrelated images is given in Table El 



on maximum likelihood and the Fisher information ma- 
trix, enabling us to estimate the expected constraints on 
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Fig. 11. The minimum variance for the parameters 
firm T, as and Us using the (Mfp)-statistics. See Fig. 
for more details. The minimum variance bound for the 
uncorrelated images is given in Table [5] 

several combinations of cosmological parameters. First, we 
assumed that only two cosmological parameters are to be 
determined from the data, fixing all the other parameters. 
In this case, using both two-point shear correlation func- 
tions ^-1- and ^_ in combination, the tightest constraints 
were obtained for the uncorrelated image configuration. 
Further, patch geometries with small cosmic variance gave 
also good results. For the aperture mass statistics (M^p), 
the best results came from a patch geometry with iV = 30 
images in 10 patches of radius R = 60'. The uncorrelated 
images could not compete with any patch geometry. 

We then took into account three and four cosmolog- 
ical parameters out of (ri,„, erg, F, $7^, 'T's)- The more pa- 
rameters are assumed to be determined from the data, 
the more important becomes large-scale information in or- 
der to resolve the near parameter degeneracies. Using the 
combined ^+ and some of the patch geometries yield 
tighter constraints than the uncorrelated image configu- 
ration. The aperture mass is best applied to patches with 

= 30 images, the results arc nearly independent of the 
patch radius in most cases. 

In most cases, the constraints obtained from the com- 
bined ^+ and ^_ were tighter than those from (M^p). 
However, the differences became smaller the more cosmo- 
logical parameters were included. 

The differences between the individual patch geome- 
tries made up to 25 percent for the minimum variance 
bound on several parameters. Thus, a 25 percent improve- 
ment on the determination on cosmological parameters 
can be obtained solely by choosing an appropriate geom- 
etry for a future cosmic shear survey. 
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